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A ESTRUTURA CELULAR DA 
CORTIÇA: LEIS ESTATÍSTICAS DE 
FORMAÇÃO DAS CÉLULAS 


M. AMARAL FORTES e M. EMÍLIA ROSA 


DEPARTAMENTO DE ENGENHARIA DE MATERIAIS, INSTITUTO SUPERIOR TÉCNICO, Av. 


Rovisco Pais, 1096 Lisboa Codex 


RESUMO 


Descreve-se a estrutura celular da cortiça, dando 
ênfase à topologia do empilhamento das células 
prismáticas. As secções da cortiça são redes 
trivalentes (três arestas em cada nó ou vértice), cujas 
arestas correspondem às intersecções com as paredes 
celulares. Analisam-se as secções do ponto de vista 
topológico e indicam-se valores experimentais das 
fracções f de células com i arestas (lados) em secções 
paralelas e perpendiculares à direcção radial. 

Em seguida procura-se modelizar a formação e 
crescimento das células na cortiça. À disposição das 
células em secções tangenciais (perpendiculares à 
direcção radial) é correlacionada com a divisão 
celular não-diferenciada do felogénio produtor das 
células. Conclui-se que apenas células com 7 ou mais 
lados sofrem divisão no felogénio, originando células- 
filhas cujas áreas diferem, no máximo, de um factor 
= 3. À modelização do crescimento na direcção radial 
permite concluir que as células de cortiça são 
geradas independentemente umas das outras, com 
alturas no intervalo 0,4h-h, em que h é a altura 
máxima das células. Aparte esta restrição, as alturas 
das células têm uma distribuição uniforme (igual 
probabilidade de todas as alturas no intervalo 
referido). 


ABSTRACT 


The cellular structure of cork is described, with 
emphasis on the topology of cell packing. The 
sections of cork are trivalent networks (three edges 
at each node or vertex), the edges of which are the 
intersections with the cell walls. The sections are 
analysed from a topological point of view and 
experimental values of the fractions f of cells with i 
edges (sides) are indicated for sections parallel and 
perpendicular to the radial direction. 

An attempt is made to model the formation and 
growth of cells in cork. The arrangement of cells in 
tangential sections (perpendicular to the radial 
direction) is correlated with the undifferentiated cel] 
division of the phellogen, the tissue that originates 
the cork cells. It is concluded that division occurs in 
this tissue exclusively in cells with 7 or more sides, 
the areas of the daughter-cells differing by a factor of 
=3, at most. The modelling of growth in the radial 
direction leads to the conclusion that cork cells are 
produced independently of each other, with heights 
in an interval 0.4h-h, where h is the maximum 
height of a cell. Except for this restriction, the cell 
heights have a uniform distribuction 
(equiprobability of cell heights in the allowed 
interval), 
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1- INTRODUÇÃO 


A cortiça pertence a uma classe de materiais 
designados por materiais celulares (ou "espumas”) 
[1], que inclui os polímeros expandidos (poliuretano 
expandido, esferovite, etc.), as estruturas em favo-de- 
mel usadas na indústria aeronáutica, e... o pão. Todos 
estes materiais são formados por células "' ocas, de 
dimensões lineares variando, conforme os casos, 
entre alguns um e alguns cm. Às células podem ser 
abertas ou fechadas, isto é, podem comunicar ou não 
umas com as outras. Toplogicamente, as células são 
poliedros, que se contactam ao longo de faces (dois em 
cada face), arestas (três em cada aresta) e vértices 
(quatro em cada vértice). Se as faces não estiverem 
completamente desenvolvidas, a estrutura é aberta. 
Se faltarem completamente as faces, a estrutura é 
formada por pequenas vigas (ou arestas), ligadas em 
nós (ou vértices). Nestes casos pode ser difícil, ou 
mesmo impossível, individualizar células 
poliédricas, mas o termo “celular” continua a ser 
aplicado a estes materiais, de que é exemplo a 
esponja natural, 


2 - ESTRUTURA CELULAR DA CORTIÇA 


Na cortiça, as células são fechadas, de forma 
prismática, e são células no sentido biológico do 
termo. São células mortas e ocas, geradas por um 
tecido - o felogénio [2] - que se encontra sob a cortiça, 
no sobreiro (Fig. la). Este tecido tem uma pequena 
espessura (uma célula, apenas) e gera células de 
cortiça para fora e células de feloderme (semelhantes 
às da madeira) para o interior. O felogénio tem 
portanto capacidade meristemática, isto é, pode 
dividir-se originando células diferenciadas (de dois 
tipos), para um e outro lado do tecido. 

Além disso, o felogénio sofre divisão celular não 
diferenciada, originando-se novas células de 
felogénio necessárias ao crescimento deste tecido. 
Quando inicialmente se forma nas extremidades do 
tronco, o felogénio desenvolve-se por divisão celular 
predominantemente segundo planos transversais 
(ver Fig. lb para nomenclatura dos planos ou 
secções), para permitir o crescimento em altura. 
Quando o crescimento em altura cessa, a divisão 
celular ocorre segundo planos radiais, para permitir o 
crescimento tangencial (isto é, o engrossamento). De 
modo que as células existentes no felogénio, numa 
dada região e num dado momento, resultaram de 
células-mães por divisão celular em duas direcções 
mutuamente perpendiculares. Estas características 
da divisão celular não-diferenciada do felogénio são 


(1 Os termos célula e celular não são usados aqui no 
sentido biológico. 


importantes se se pretender modelizar o crescimento 
deste tecido, o que faremos neste artigo. 

As células de cortiça são empilhadas em colunas, 
de eixo paralelo à direcção radial da árvore (Fig. 1b). 
Cada célula do felogénio gera, uma a uma, as células 
de cortiça de uma coluna. 

Numa secção perpendicular à direcção radial - 
designada por secção tangencial - observa-se uma 
disposição do tipo "favo-de-mel", com células 
poligonais, relativamente regulares (Figs. lb e 2a). 
Esta disposição é idêntica à das células do felogénio 
visto que cada célula deste tecido origina uma coluna 
de células de cortiça. O número de lados de cada 
célula poligonal é variável (Fig. 2a). Mas em cada 
vértice encontram-se sempre três células e, portanto, 
três arestas (lados das células), o que significa que a 
rede bidimensional observada na secção tangencial é 
uma rede trivalente (três arestas em cada vértice). 


na SECÇÃO 
rs TANGENCIAL 
cortiça Felogénio 4 q JANGENC TRANSVERSAL 
Es - “oADIAL se: Ss” aos Era 
| Perto Sings rg 0 00h 
| Marli É H, LIA 
| | URL 
| | UPA Tn TA 
a E po DE Aa AP 
L SECÇÃO INCC A 
MS RADIAL “ITA” | SECÇÃO 
Nº | TANGENCIAL 
(a) (b) 


(c) 


FIGURA 1 


Fig. 1-a) Esquema do tronco do sobreiro, 
mostrando a camada exterior de cortiça, gerada pelo 
felogénio. b) Esquema da estrutura celular da cortiça 
mostrando as células prismáticas (não 
necessariamente hexagonais) empilhadas em 
colunas, com bases desencontradas. Indica-se 
também a nomenclatura usada para as direcções e 
para as secções. c) Uma célula é contactada nas faces 
laterais pelas bases das células vizinhas; dimensões 
típicas: H=35um,a =14ym. 


As células em colunas adjacentes (geradas por 
células adjacentes do felogénio) não têm as bases no 
mesmo plano (Figs. 1b, c e 2b). Em princípio, cada 
célula do felogénio gera células de cortiça de uma 
forma independente da das células de felogénio 
adjacentes, de modo que não há correlação entre as 
alturas das células de cortiça adjacentes. O resultado 


4 TÉCNICA 2 - 88 


(a) 100 um (b) 
FIGURA 2 
Fig. 2 Secções da cortiça observadas no 


microscópio electrónico de varrimento: a) secção 
tangencial; b) secção radial. 


é um empilhamento desencontrado nas várias 
colunas (Figs. 1b, c), como se pode verificar na secção 
radial da Fig. 2b. Esta secção contém a direcção 
vertical, ou axial (eixo da árvore), e a direcção radial. 
Mas qualquer secção que tenha a direcção radial tem 
o mesmo aspecto. À direcção radial é uma direcção de 
simetria da estrutura da cortiça. Designaremos por 
secções radiais todas as que contêm a direcção radial. 
Todas estas secções têm a estrutura em ladrilho que 
se pode observar na Fig. 2b, e que pode ser descrita 
(embora com alguma “agressão” à topologia) como 
um empilhamento desencontrado de rectangulos. À 
“agressão” à topologia tem a ver com o facto de, 
nestas secções, cada célula ter um número de lados (e 
vértices) geralmente superior a 4, não sendo portanto 
rectangular no sentido topológico. De facto, tal como 
na secção tangencial, cada célula é um polígono com 
um número de lados e vértices variável (Fig. 2b), mas 
há sempre três arestas em cada vértice. Quer dizer, 
topológicamente as duas secções são muito parecidas, 
embora geometricamente seja notória a diferença 
entre elas. 

Um aspecto importante da estrutura celular da 
cortiça é o ondulamento das paredes laterais das 
células, observável na Fig. 2b [3, 4). As bases das 
células também não são planas, mas em geral não 
apresentam ondulações. Estas irregularidades são 
resultado do modo de crescimento da cortiça, com as 
novas células geradas pelo felogénio a ter de 
empurrar as células já formadas. As paredes 
celulares sofrem flexão e encurvadura (mas não 
fractura!) devido a estas tensões de crescimento. Mas 
a verdade é que as células são geradas com o 
tamanho correcto para poderem ser direitas e planas. 
Quer dizer, se não fosse a necessidade de empurrar 
células já formadas, a cortiça teria células de paredes 
planas, paralelas e perpendiculares à direcção radial. 
Que isto é assim, infere-se da Fig. 3, onde se mostram 
as secções tangencial e radial de cortiça após 
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aquecimento durante lh a 300º C, ao ar. O 
aquecimento provoca o endireitamento das paredes, 
ao mesmo tempo que a cortiça expande e se removem 
as tensões internas de crescimento. 

Para completar esta breve descrição da estrutura 
celular da cortiça, convém referir ainda que as 
dimensões das células variam com a época do ano em 
que se formam, mas são, em média, as seguintes [4]: 
altura 30-40 um; aresta de base: 13-15 um; espessura 
da parede celular: 1-2um. As células de cortiça são 
mais pequenas do que as dos materiais celulares 
artificiais. 

As células poliédricas (prismáticas) na cortiça 
contactam as células adjacentes através das faces, 
sendo o número médio de células adjacentes de uma 
dada célula =14, como se pode inferir da descrição 
feita (Fig. 1c). (De facto o número 14 é o resultado 
correcto se não houver qualquer correlação entre a 
altura e o número de faces laterais das células [4]). 

Neste artigo procuraremos modelizar a divisão 
celular do felogénio de modo a que conduza às 
distribuições f (fracção de células com i lados) 
observadas nas secções tangenciais e radiais da 
cortiça (Tabela 1). As "leis" de divisão celular que 
terão de ser introduzidas na modelização para obter 
as distribuições correctas, são, como veremos, leis 
plausíveis, mas não será feita qualquer tentativa de 
explicação física ou biológica da necessidade dessas 
mesmas leis. Não obstante, a análise que iremos 
fazer permite descortinar algumas características 
“médias” dos fenómenos de divisão celular não- 
diferenciada e diferenciada, as quais são de utilidade 
para uma análise mais aprofundada desses 
fenómenos. 

Preocupar-nos-emos exclusivamente com a 
topologia das células e não com a sua geometria; em 
particular, as ondulações das paredes serão 
irrelevantes na discussão que se segue. As 


(a) 100 jum (b) 


FIGURA 3 
Fig. 3 Efeito do aquecimento na estrutura da 


cortiça. Compare-se com a Fig. 2 e note-se o 
endireitamento das paredes celulares: a) secção 
tangencial; b) secção radial, 


características geométricas são igualmente 
importantes e poderiam ser estudadas de uma forma 
semelhante ao que se fará para as características 
topológicas. São no entanto mais difíceis de medir, 
particularmente devido à irregularidade geométrica 
das células, resultante das ondulações das paredes, 


3 - TOPOLOGIA DAS SECÇÕES 


Em qualquer secção as células são poligonais, 
com um número de lados (arestas), ou de vértices, 
variável de célula para célula. Nas secções radiais, 
em ladrilho, o número de vértices (e portanto de 
arestas) de uma célula inclui os vértices em que se 
juntam a essa célula as arestas das células em 
colunas adjacentes. Na Fig. 4 indicam-se os números 
de lados (ou vértices) de células em redes que são 
geometrica e topologicamente semelhantes às que se 
observam nas secções tangenciais e radiais. 


FIGURA 4 


Fig. 4 Redes trivalentes, análogas às que se 
observam nas secções da cortiça: a) tangencial;, b) 
transversal ou radial. Está indicada a poligonalidade 
de cada célula. 


Como já se referiu, em qualquer secção os 
vértices são trivalentes, isto é, há 3 arestas em cada 
vértice. Por isso, se for À o número de arestase Vo 
número de vértices, tem-se: 


24=3V (1) 
Por outro lado, a equação de Euler relaciona o 
número F de células (ou faces, na designação 


topológica) com A e V: 


F+V=A (2) 


desde que o número de células seja grande. 
Combinando as duas equações anteriores obtém-se 
2F = V; como há 3 células em cada vértice, o número 
médio de vértices por célula é 3V/F=6. Isto é, 
tomando um grande número de células e 
determinando a fracção f de células com i lados (i- 
células), tem-se: 


2.151 
(3) 
<i>= > if,=6 


Uma estrutura celular bidimensional com 
vértices trivalentes terá a máxima ordem topológica 
se todas as células tiverem 6 lados. Na Fig. 5a,b 
indicam-se duas variantes geométricas desta 
estrutura com f=1, que são, respectivamente, 
semelhantes à secção tangencial e às secções radiais. 
As redes da Fig. 5a,b têm uma considerável ordem 
geométrica. Na Fig. 5c,d mostram-se estruturas 
semelhantes , ainda com a máxima ordem 
topológica mas sem qualquer ordem geométrica. 


FIGURA 5 


Fig. 5 Redes trivalentes com todas as células 
hexagonais: a) e b) têm ordem elevada; c) e d) são 
desordenadas. 


A desordem topológica de uma rede bidimensional 
trivalente pode ser medida pelo segundo momento, 
ua, da distribuição f/, em relação à média: 
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p= > (1-6)/f = (1-6) 
l 
(4) 


A ordem máxima corresponde a u3=0. Em redes 
bidimensionais ny pode variar entre zero e infinito. 


TABELA 1 
Distribuição de poligonalidades nas secções da 
cortiça 


O | rensencal | mada o 


Na Tabela 1 indicam-se as distribuições f; para os 
dois tipos de secções, e os respectivos uy. Estes 
valores foram obtidos numa única amostra, mas 
medidas noutras amostras indicam que eles são 
característicos da cortiça. Os valores relativamente 
baixos de 42 reflectem uma certa ordem no arranjo 
das células, que iremos estudar com mais pormenor. 
E de notar a ausência de células com 2 ou 3 lados e a 
raridade de células com 9 ou mais lados. É de referir 
ainda que as células com mais lados são em geral as 
de maior dimensão, mas não sistematicamente. 


4 - SECÇÃO TANGENCIAL: DIVISÃO NÃO- 
DIFERENCIADA DO FELOGÊNIO 


Como se referiu já, a disposição das células de 
cortiça na secção tangencial reflecte a disposição das 
células do felogénio que as produz. O felogénio é um 
tecido bidimensional, que cresce por divisão celular 
(não-diferenciada), isto é, por divisão das células "ao 
meio”, acompanhada de crescimento celular. Esta 
divisão nada tem a ver com a divisão diferenciada 
que gera as células de cortiça. Do ponto de vista 
topológico, há que procurar um modelo de divisão 
celular que conduza a um tecido com a distribuição f 
da secção tangencial da cortiça, indicada na Tabela 
1. Uma célula com i lados origina por divisão duas 
células com i e i, lados, sendo (Fig. 6): 


i+i,=i+4 (5) 
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Pode concluir-se que no felogénio não ocorrem 
divisões conducentes a células com 3 ou 2 lados (Fig. 
6b,c) visto que as células deste tipo não são 
observadas na secção tangencial da cortiça. 


(a) (b) (c) (d) 
FIGURA 6 
Fig. 6 Exemplos da divisão celular em células 


bidimensionais. O número de lados i;, io, das células- 
filhas está indicado, bem como o número i de lados da 
célula-mãe (número com círculo). É sempre i=ij +ig- 
4. No felogénio não aparecem células com dois ou três 
lados (respectivamente b e c), nem há divisão de 
células com menos de 7 lados. 


Por outro lado, devem ser raras as divisões de células 
com 4 lados, as quais conduziriam a duas células 
adjacentes com 4 lados (Fig. 6d), o que também não se 
observa na secção tangencial. 

Na modelização, ou simulação, que iremos fazer 
da divisão celular não-diferenciada do felogénio, 
admitiremos primeiro que a divisão ocorre apenas 
em duas direcções mutuamente perpendiculares 
daquele tecido, o que se coaduna com as 
características da divisão celular mencionadas na 
Introdução. Consideremos então um quadrado de la- 
-do unitário (Fig. 7) e tomemos nele um ponto P ao 
acaso, escolhendo um par de números aleatórios x, y, 
no intervalo (0, 1), para coordenadas do ponto (todos 
os valores no intervalo (0, 1) são equiprováveis). O 
quadrado é dividido em dois por um segmento 
passando por P e paralelo a x ou y, ao acaso (igual 
probabilidade de uma e de outra): Fig. 7a. 

O processo de fragmentação continua da mesma 
forma (Fig. 7b), tomando sucessivos pontos ao acaso € 
dividindo, paralelamente a x ou y, a célula onde o 
ponto se localiza. Por este processo atingir-se-á 
eventualmente um estado estacionário, isto é, uma 
situação em que não há alteração expectável de f; 
quando ocorrem divisões de mais células. Na Fig. Tc 
mostra-se o resultado de uma fragmentação de 
acordo com este modelo, e contendo 50 células. À 
distribuição f, correspondente (Tabela 2, coluna B) 
está muito longe da que se observa na secção 
tangencial da cortiça (Tabela 2, coluna A). De facto, e 
se bem que não haja células com menos de 4 lados (o 
que é inerente ao modelo), a fracção de 4- células é 
demasiado grande, como se pode ver na Tabela 2, 
comparando as colunas À e B. Neste modelo a 


FIGURA 7 


Fig. 7 Divisão de um quadrado por sucessivas 
divisões celulares: a), b), indicam o modo de efectuar 
as divisões, após localização de um ponto P, ao acaso; 
c) mostra o resultado de 50 divisões, sem qualquer 
restrição quanto às células que sofrem divisão e ao 
modo como se dividem; d) resultado de = 200 divisões 
com as restrições indicadas no texto, o que conduz a 
uma distribuição de poligonalidades muito 
semelhante à da secção tangencial da cortiça. 


probabilidade de divisão celular é proporcional à área 
das células, o que é plausível, mas a divisão ocorre 
com equiprobabilidade dos tamanhos das duas 
células-filhas, o que é pouco plausível. Na verdade, é 
de esperar uma divisão mais ou menos a meio. 

Este modelo de divisão celular pode ser 
melhorado, introduzindo uma ou mais das seguintes 
restrições: 

TABELA à? 


Distribuição de poligonalidades na secção tangencial, obtidas com diferentes masdelos de divenho eslular 


Identificação das colunas simalações do tipo da Fig. 7-d; Simulação do tipo da Pig. Ea 


dc Valores experimentais C -re6;f= 10; ne 150 FarT;f=18,n=1D0 
Seumulação do tipo da Fig, T-e D-de7,f=10,n=150 Simulação do tipo da Pig. É-b 
B Divisão celularemo ed;f=], nf de cólolaa, n= Bi) E -set;f=17;n=>200 Gr ETi=1820=100 


1 - Tamanho mínimo para que ocorra divisão; este 
tamanho tem de ser definido em relação à área média 
das células, uma vez que no modelo, o crescimento 
das células está “congelado”. 


2 - Restrições relacionadas com o número de lados das 
células que podem dividir-se, p. ex., admitir que só se 
dividem células com 6 ou mais lados. 

3 - Divisão segundo a menor dimensão, isto é, o 
segmento introduzido na divisão é paralelo à menor 
dimensão da célula. 


4 - O segmento que origina a divisão situa-se na 
região central da célula, cuja largura é uma fracção, 
f, da largura total da célula. Admitiremos que todas 
as divisões nesta região são equiprováveis. As 
células-filhas têm áreas que diferem, quanto muito, 
de um factor (1+1)/(1-f). 


Foram efectuadas simulações da divisão celular, 
introduzindo várias restrições do tipo acima 
indicado. Em todos os casos foram geradas 100 
células e determinados os f,. Nalguns casos 
continuou-se a simulação até 150 células, não se 
verificando variações significativas dos f,. Os valores 
apresentados podem pois considerar-se 
representativos de um estado estacionário. 

Na Tabela 2 indicam-se alguns exemplos de 
resultados obtidos em diversas condições, indicadas 
junto da Tabela. Verifica-se que os valores da coluna 
E são bastante próximos dos valores observados na 
secção tangencial da cortiça (Tabela 2, coluna A): 
Aqueles valores foram obtidos nas seguintes 
condições: 


a) escolha da célula que sofre divisão - probabilidade 
proporcional à área, isto é, toma-se um ponto ao 
acaso e efectua-se a divisão na célula em que o ponto 


se situa; 
b) divisão só em células com 7 ou mais lados; 


c) divisão perpendicularmente à maior dimensão, 
ocorrendo na região central com f= 1/2. 


Na Fig. 7d mostra-se o resultado de uma 
simulação nestas condições. A distribuição f 
correspondente figura na coluna E da Tabela 2. 

Admitiu-se até agora que a divisão celular pode 
ocorrer em duas direcções mutuamente 
perpendiculares, com igual probabilidade. Para ver 
até que ponto esta restrição influencia os valores 
finais de f, efectuaram-se duas outras simulações em 
que apenas se alterou a condição relativa às direcções 
de divisão celular (isto é, as outras condições, a, b, c, 
acima indicadas foram mantidas; a maior dimensão 
de uma célula é a maior corda da célula). Na primeira 
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(a) (b) 


FIGURA 8 


Fig.8 a)Divisão celular segundo a direcção 
perpendicular à maior corda da célula que se divide; 
b) divisão celular exclusivamente numa direcção. 


destas simulações fizeram-se sucessivas divisões dum 
círculo , o que conduz a uma distribuição aleatória da 
orientação das paredes celulares (Fig. 8a). Na outra 
simulação, admitiu-se apenas uma orientação 
possível de divisão, tendo-se partido de uma rede de 
rectas perpendiculares à direcção de divisão (Fig. 8b). 
Em qualquer das duas simulações as distribuições fi; 
resultantes ao fim de 50-100 divisões são muito 
próximas das obtidas com divisão em duas direcções 
perpendiculares. As distribuições f; correspondentes 
figuram nas colunas F e G da Tabela 2, 
respectivamente. Conclui-se que a topologia 
resultante da divisão celular não é afectada pela 
geometria inicial nem pela orientação em que é feita 
a divisão. Em particular, a frequência (ou peso) 
relativa das divisões celulares no felogénio, na 
direcção axial e na direcção tangencial, não tem 
incidência sobre a topologia resultante. 


5 - SECÇÕES RADIAIS 


À estatística topológica do arranjo celular nestas 
secções é mais fácil de modelizar. Cada célula do 
felogénio produz células de cortiça, uma após outra. 
Admitamos que as células têm alturas no intervalo 
(£, 1)* e que todas as alturas neste intervalo são 
equiprováveis. As células vão sendo produzidas com 
alturas ao acaso. Esta hipótese é razoável, embora 
seja possível a existência de uma correlação entre a 
altura média das células de uma coluna e a área da 
célula de felogénio que as origina. Admitiremos, pois, 
que a probabilidade p(x)d(x) de uma célula de cortiça 


*. Isto é, o cociente entre a altura mínima e a altura 
máxima é E. A distribuição de poligonalidades nas 
secções paralelas à direcção radial depende então 
apenas de £ e não do valor da altura máxima a qual 
pode considerar-se unitária 


ter uma altura no intervalo x, x+ dx é uniforme no 
intervalo [6,1], isto é: 


plx) =0 q<Ex>1 


O problema do arranjo topológico das células reduz- 
se a um problema de probabilidade geométrica. Há 
uma distribuição D de segmentos adjacentes numa 
recta (Fig. 9a), com uma distribuição p(x) de 
comprimentos (x< 1). 


a ? (a) 


(Db) 


pese fada e 
po Le 
ame (c) 
pe | —el 
FIGURA 9 


Fig.9 a) Distribuição de segmentos numa recta, 
com uma distribuição p(x) de comprimentos x. L é um 
segmento colocado sobre a recta, contendo 6 ponto-E, 
b) Uma célula da fiada central tem um número de 
lados que depende do número (ny ng ) de pontos E, em 
fiadas adjacentes, dentro da célula. .c) Definição de 
L* e de outras quantidades usadas na dedução do 
Apêndice. 


As extremidades destes segmentos são designadas 
por pontos-E. Consideremos um segmento de 
comprimento L(L<=1) colocado ao acaso na recta 
(todas as posições da extremidade esquerda, por 
exemplo, de L são equiprováveis). Qual é a 
probabilidade de L conter exactamente n pontos-E? 
Designaremos essa probabilidade por P(n, L). No 
apêndice calcula-se P(n, L), relacionando-a com p(x), 
egs. 21-23. 


Voltando ao problema topológico das secções 
radiais, consideremos uma célula de altura L que 
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contém (n,, n,) vértices adicionais de cada lado (Fig. 
9b), resultantes do contacto com células adjacentes. 
A poligonalidade i dessa célula é: 


(7) 


O número expectável de células contendo (ny, no) 
vértices adicionais é 


l 
q 
E: | o, 1) P(n plJjdL 


(8) 


coma=1paran én,ea=2paran =n,. À fracção de 
células com i lados obtém-se finalmente somando as 
quantidades (8) para todas as decomposições (n,, n,) 
de i compatíveis com a eg. 7. 

No caso da cortiça verifica-se que não há células 
com mais de 10 lados, o que significa não haver mais 
do que 3 vértices adicionais de cada lado, ou seja, 
€> 1/3. Por outro lado, a existência de células com 8 
lados implica £<1/2. Os cálculos foram feitos para 
valores de € neste intervalo, indicando-se na Tabela 3 
os resultados obtidos. Para £= 0,4 a concordância com 
os valores experimentais (primeira coluna da Tabela) 
é excelente. 

Conclui-se que as células têm alturas 
uniformemente distribuidas no intervalo 0,4-1: o 
quociente entre a altura mínima e a altura máxima 
das células é portanto 0,4. Parece também poder 
concluir-se que não há correlação entre as alturas das 
células adjacentes numa mesma coluna, nem entre as 
alturas de células adjacentes em colunas adjacentes, 
isto é, a altura das células não é correlacionável com 
o tamanho ou poligonalidade das células de felogénio 
que as produzem. 


TABELA 3 


Valores experimentar e calculados de f, na secção radial 


Valores 
1! ss 
esperimentaim 


0,0367 0, 0288 


OD. 2000) 0,1857 
0,5271 | 05H11 


0,1992 D,1856 


0,0366 


0,0268 
0.0003 19105 


13u 10 1,3n10* 


6 - CONCLUSÕES 


A cortiça é formada por células 
(geometricamente) prismáticas, empilhadas em 
colunas paralelas à direcção radial Cada célula deve 
ser encarada, do ponto de vista topológico, como um 
poliedro com faces, arestas e vértices. O número 
médio de faces de cada poliedro (igual ao número 
médio de células adjacentes a uma dada célula) é 
= 14, Em cada aresta há três faces e em cada vértice 
há quatro arestas. À estrutura tridimensional é 
tetravalente (quatro arestas em cada vértice). Nas 
secções da cortiça, as células são poligonais com um 
número médio de lados igual a 6, e com três arestas 
em cada vértice (redes trivalentes). Embora todas as 
secções tenham esta topologia, elas são, do ponto de 
vista geométrico, diferentes umas das outras. Na 
secção tangencial (perpendicular à direcção radial) a 
estrutura é semelhante a um favo-de-mel, embora 
menos regular. Nas secções contendo a direcção 
radial (secções radiais), as células têm uma forma 
(geométrica) rectangular, e a estrutura lembra um 
ladrilho. Em qualquer dos casos, o número de lados 
das células (cujo valor médio é igual a 6) é variável de 
célula para célula, mas as fracções f de células com i 
lados são características de cada secção, variando 
pouco de amostra para amostra. Esta constância tem 
certamente a ver com a existência de "regras rígidas" 
relacionadas com o modo como são geradas as células 
na cortiça. 

À distribuição f em secções tangenciais pode ser 
identificada com a distribuição f, das células de 
felogénio que as originam. Por sua vez, estas formam- 
se por divisão celular segundo planos radiais, sendo 
possível estabelecer um modelo de divisão que 
conduz à distribuição f; experimentalmente 
observada. Este modelo caracteriza-se por divisões 
que originam células-filhas com áreas semelhantes 
(diferindo no máximo de um factor 3) e por ausência 
de divisão em células com menos de 7 lados. A 
probabilidade de divisão celular é , âparte esta 
última restrição, proporcional à área das células. 

Por outro lado, a distribuição f; em secções 
contendo a direcção radial pode ser explicada 
admitindo que as células de cortiça produzidas pelo 
felogénio têm uma distribuição uniforme de alturas 
no intervalo 0,4h-h, em que h é a altura máxima das 
células. As células são produzidas com alturas ao 
acaso dentro deste intervalo, sem correlação entre as 
alturas de duas células produzidas em sucessão ou 
entre as alturas de células adjacentes em diferentes 
colunas. 

Em conclusão, a distribuição aparentemente 
desordenada das células nas secções da cortiça 
(diferentes tamanhos, diferentes números de lados, 
paredes (arestas) curvas) esconde as leis de formação 
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das células, as quais podem ser deduzidas por uma 
análise cuidada e são, do ponto de vista estatístico, 
leis relativamente simples. Particularmente 
importante é a possibilidade de modelizar a divisão 
celular não-diferenciada, deduzindo leis 
probabilísticas relacionadas com a topologia das 
células. 

Fica infelizmente sem resposta, no âmbito desta 
análise, a seguinte questão fundamental: porque é 
que as células com mais lados são as mais adultas, 
no sentido de mais reprodutivas; porque é que estas 
se tornam instáveis e criam uma parede celular que 
as divide em duas, ao mesmo tempo que o núcleo 
sofre mitose para dar os núcleos das duas células- 
filhas. Em suma, quais são as "forças" responsáveis 
pela divisão celular? 


APÊNDICE 


Distribuição de poligonalidades em secções radiais 


Consideremos uma recta dividida em segmentos, 
com uma distribuição de comprimentos, x, definida 
por p(x), tal que p(x)dx é a probabilidade de um 
segmento ter comprimento entre x e x+dx (Fig. 9a). 
Não há qualquer correlação entre os comprimentos 
de segmentos adjacentes. Os comprimentos podem 
variar no intervalo (£, 1). O caso em que o intervalo 
permitido é infinito, como acontece com uma 
distribuição de Poisson de pontos sobre uma recta, 
pode ser tratado da mesma maneira e será abordado 
no fim deste Apêndice. Os pontos-extremidade dos 
segmentos serão designados por pontos-E. Tem-se: 


1 
l =| xp(x)dx 
á 


l 
| plx)dx=1 : 
Ê (9) 


em que [é o comprimento médio dos segmentos. 

Determinemos em primeiro lugar a 
probabilidade m,(L)dL de que a soma dos 
comprimentos de n segmentos (nz 1) esteja no 
intervalo L, L + dL. Consideraremos apenas 


ESL<I (10) 
Obviamente 


n(L)=p(L) (1) 


Para n=2 (dois segmentos), um deles tem 
comprimento x e o outro L-x, com qualque x no 
intervalo £, L-&. Tem-se 


L—t 
ndlo=| plo) p(L — x)dx 2€<=L<l 
E 
(12) 


n(L)=0 L<2€ 


Usando o mesmo raciocínio é possível obter a 
seguinte fórmula de recorrência 


L=nt 
Mo pqt)= | p(x)n, (L — x)dx (n+1I)jf<SL<i 
2 


(13) 


14 (1)=0 L<(n+1)% 

A probabilidade de que a soma dos comprimentos de 
n segmentos seja inferior a L(L<1), que 
designaremos por p (L), é: 


L 
p,L)= | n,(x)dx né<sL<] 
né, 
(14) 
pXL)=0 L<nt 
sendo 
p(L)=1 


Consideremos agora um segmento de 
comprimento L, e determinemos a probabilidade 
P(n,L) de que existam, em L, exactamente n pontos- 
E. Admitiremos, como atrás, que E<L<1. A 
extremidade esquerda do segmento L pode 
encontrar-se num segmento de comprimento [, l+dli, 
com probabilidade proporcional a lp(l)dl. A 
probabilidade é pois (1 Dp(Ddl. 

Para que haja O pontos em L é preciso que a 
extremidade esquerda de L esteja num segmento 
[>L, numa região de comprimento [-L (de modo a 
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que a extremidade direita se situe também dentro de 
[), ao que corresponde a probabilidade (l-L)/l. Obtém- 
se então 


l 
P(O,L)= | u-L)PP a 
lt 1 


(15) 


Para calcular a probabilidade P(n,L) com n>o0, 
definamos o comprimento de excesso L*, que é o 
comprimento do segmento JL que ultrapassa a 
extremidade do segmento [ onde L tem a sua 
extremidade esquerda (Fig. 9c): 


L*=L+yl (16) 


onde y designa a distância entre as extremidades 
esquerdas de L e |. Para que L contenha exactamente 
n pontos (n>0) é necessário que L>(n - 1) €e que L* 
exceda a soma dos (n-1) segmentos que se seguem a Í, 
mas não exceda a soma dos n segmentos que se 
seguem a [. De modo que 


 p(D 
Penta | | Cal|p, (Lty-D-p,d:+y- Do 
Ly 


(17) 


Se | 2L os limites de y são LL e |; se I<I, esses 
limites são 0,/. Por outro lado p (L*) anula-se se 
L*<kt. Pode concluir-se que os intervalos que 
contribuem para o integral da eq. 17 são os que estão 
indicados na Tabela 4; estes intervalos dependem do 
valor de L, sendo necessário considerar três 
intervalos de L. 


TARELA 4 


Limites de integração em Lo y na eg. 17 


=== 


“Significa que não há nenhum intervalo 


E Lonk 


L-nk : L-fn-1)E 


LAn-1H; E 


Lj;l 


a=[-L+nk 
b=L+(a-1E 


Intervalo O ma 
de ; 


É agora possível efectuar os cálculos conducentes 
a P(n,L) para uma qualquer distribuição p(x) dos 
comprimentos dos segmentos. Para uma distribuição 
uniforme dos comprimentos no intervalo €,1 (cf.eq.6): 


1 
o A Esy<l 
(18) 
pl)=0 EE + gl 
obtém-se: 
(L- ng") 
(LL) = (= 
F (n-19)!(1- 
Da 19) 
(L—n&” 
p (L)= ——— nE<L<l 
i n1 8" 
(20) 


sendo nulos os valores fora do intervalo de L indicado. 
A eq. 20 é válida para qualquer n. Utilizando estas 
expressões e a Tabela 4 é possível calcular P(n,L). 
Para P(0,L) obtém-se da eq.15: 


|=15 
pro L)= 
me 2 
(21) 
Definamos 
(L — ng” 
“In IJ =" 
AO em 


Obtém-se, para n>1,a partir da eq. 17: 


2 
1=L2(n+1E P(nl)=-—(Q, ,4-20,+0.,) 
E=É == ! n 


-20,) 


n-1 n 


(n+1EZL=>nE Pl 6 
1 E 


(23) 


nE=L>(n- DE Pln,l)=Q, 


12 TÉCNICA 2 - 88 


O problema da distribuição de poligonalidades 
em secções radiais pode ser assim enunciado: é dada 
uma distribuição p(x) de segmentos e duas outras 
distribuições p(x) não correlacionadas com a 
primeira (Fig. 9b). Qual é a fracção expectável de 
segmentos da primeira distribuição que contêm n, 
pontos da segunda distribuição e n, pontos da 
terceira distribuição? Essa fracção, (n n,), é 


l 
(n Po)= |, p(L) P(n ,L)P(n,,L)dL 


(24) 


donde resulta a eg. 8 do texto. 

Para finalizar este Apêndice consideremos uma 
distribuição de Poisson de pontos sobre uma recta. À 
distribuição dos comprimentos x dos segmentos 
(O<x<o) é [6] 


= a a 
x plx)=e (25) 


e a probabilidade de n pontos em L(L>0) é, para 
qualquer n>0, 


| (Ux e Ria 
Pí(inL) = 
n! 


(26) 


que é um resultado bem conhecido [6]. A fracção de 
células com n,, n, pontos adicionais é dada pela eq. 24 
com limites de integração 0,. Obtém-se 


| 
(n, + no)! 
(n ny= 
l nj+no+1 
3 n no! 
(27) 
A fracção de células com i lados obtém-se somando 
estas quantidades para todos os valores compatíveis 
com a eq. 7, com o seguinte resultado: 


(28) 


que verifica as identidades (3) do texto. Os valores de 
j, obtidos de (28) estão indicados na última coluna da 
Tabela 4. 
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RESUMO 


A Modelação pode ser definida como a 
descrição de um processo tecnológico, partindo dos 
fenômenos físicos em que se baseia, e usando uma 
linguagem matemática para o descrever 
quantitativamente, por forma a poder optimizá-lo. 
Após introduzir o conceito de modelação (31), indica- 
se o papel que desempenha nas várias fases de um 
processo de desenvolvimento tecnológico (82), o que 
permite evidenciar a relação entre as ciencias 
fundamentais (matemática e física) e a sua aplicação 
ao estudo de obras de engenharia (83). A partir desta 
é feita uma lista de tópicos relevantes (Tabela 1) para 
uma Área Cientifica de Modelação (84), conduzindo a 
um elenco de disciplinas (Tabela II), cuja inserção nos 
planos de estudo vigentes no 1.S.T. é delineada (85). A 
experiência adquirida, e a actividade em curso, no 
âmbito de um dos departamentos do 1.S.T. (86), serve 
de base para algumas conclusões (87) sobre as 
possibilidades de desenvolvimento futuro. À 
actividade de ensino de modelação é testemunhada 
pela publicação de 2 livros, 1 por entidade nacional e 
1 por editora internacional, que estão incluidos 
(referências 1-2) na bibliografia breve, anexa ao 
presente artigo. 
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$ 1 INTRODUÇÃO 


A Tecnologia é a aplicação da ciência à 
realização de tarefas práticas, e como tal, as obras de 
engenharia são uma extensão, concebida pelo 
homem, dos processos que ocorrem na natureza. Da 
definição de tecnologia que se acaba de apresentar, 
decorre o papel fulcral desempenhado pela 
modelação, como o elo que permite aplicar os métodos 
científicos à resolução de problemas de engenharia, 
por exemplo, para verificar se um processo proposto é 
praticável, ou optimizar o projecto de um 
equipamento. A modelação é também essencial para 
a ciência, como a descrição quantitativa de 
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fenómenos que ocorrem na natureza; a tecnologia 
alarga substancialmente o âmbito da modelação, pelo 
estudo de processos artificiais, concebidos pelo 
homem, para modificar a natureza ou superar as suas 
limitações. À medida que o progresso tecnológico se 
acelera, e a sua ligação com a ciência e a engenharia 
se aprofunda, a modelação adquire uma importância 
sempre crescente, como forma de viabilizar a 
realização de equipamentos sofisticados e 
concretização de obras complexas. A tecnologia é o 
simbolo mais saliente do século XX, isto é, o dominio 
em que a sociedade contemporânea mais progrediu, e 
melhor se distingue do passado; a modelação, como 
componente intelectual da tecnologia, é porventura o 
ramo do conhecimento humano que mais evoluiu no 
presente século, e melhor caracteriza a cultura 
contemporânea. 


$ 2-PROCESSO DE DESENVOLVIMENTO 
TECNOLOGICO 


Para documentar o papel da modelação 
matemática e física no progresso tecnológico, 
indicamos a sua intervenção nas várias fases de uma 
obra de engenharia: (i) o desenvolvimento de um 
equipamento visa satisfazer a um conjunto de 
especificações técnicas, que resultam de um estudo 
das necessidades do mercado, isto é, de um modelo de 
utilização dos equipamentos actuais, que permite 
prever as vantagens de introduzir certos 
melhoramentos na geração seguinte; (ii) o projecto 
preliminar baseia-se num modelo global do 
equipamento, envolvendo parâmetros que afectam a 
sua configuração, capacidade, complexidade e custos, 
e devem ser escolhidos de forma a maximizar a 
eficiência global; (iii) o projecto detalhado envolve 
modelos mais pormenorizados de cada sub-sistema, 
usados para garantir que este dá o contributo 
esperado para o equipamento global, e para estudar 
problemas de concretização de cada parte; (iv) a 
construção e experimentação de protótipos é 
precedida e acompanhada por simulações usando 
modelos, por forma a detectar e explicar 
discrepâncias, antes de estas se transformarem em 
acidentes; (v) a certificação e validação do 
equipamento consiste de um conjunto de ensaios, que 
podem ser conduzidos sem risco a ritmo acelerado, 
desde que o resultado em cada ponto de ensaio 
coincida com a previsão do modelo dentro de certas 
tolerâncias; (vi) estando o modelo validado por 
comparação com a experiência, pode ser usado para 
simular condições de funcionamento marginais, que 
seria perigoso concretizar, ou para testar um grande 
número de configurações, das quais só algumas 
podem ser concretizadas, por razões de eficiência e 
custo; (vii) o planeamento da produção obedece a um 


modelo de aprovisionamento e optimização do 
processo, e a fabricação automatizada com elevada 
precisão depende de um modelo detalhado de cada 
componente; (viii) a operação de um equipamento 
complexo e sofisticado é em geral acompanhada por 
um simulador, que serve para treinar os operadores, 
experimentar condições anormais de funcionamento 
e antever as vantagens de certos melhoramentos; (ix) 
a manutenção do equipamento utiliza indicadores de 
usura, que permitem substituir apenas os 
componentes gastos, antes de se verificarem avarias, 
e sem desperdiçar recursos com revisões 
desnecessárias; (x) o modelo do equipamento pode ser 
utilizado em caso de avaria, para localizar 
rapidamente o componente deficiente, evitando que 
se paralise o equipamento seguindo pistas falsas; (x1) 
o aprovisionamento de peças de reserva e 
planeamento do trabalho de manutenção, recorre a 
modelos de usura, que permitem manter o 
equipamento em serviço, sem empatar capital em 
infraestruturas de apoio sobredimensionadas. 


83 - RELAÇÃO ENTRE FUNDAMENTOS 
E APLICAÇÕES 


A modelação, para ser quantitativa, baseia-se 
em ciências fundamentais, como a matemática e a 
física, mas para que seja relevante a uma 
especialidade de engenharia, deve ser dirigida por 
um conhecimento aprofundado desta. Assim, os 
quatro contributos essenciais para a modelação, 
correspondem a: (1) do ponto de vista matemático, o 
cálculo exacto, aproximado ou optimizado, em termos 
de variáveis que descrevem o funcionamento de um 
certo equipamento; (ii) do ponto de vista físico, a 
representação realista dos fenómenos ou processos 
em que se baseia o funcionamento do equipamento, 
incluindo os efeitos importantes e omitindo tudo o 
que seja acessório; (iii) a selecção de métodos 
matemáticos e lista de efeitos físicos a incluir 
baseiam-se no conhecimento da especialidade, por 
exemplo, das condições em que o equipamento vai ser 
utilizado; (iv) a interpretação dos resultados do 
modelo depende dos objectivos e desempenho que se 
pretende do equipamento, e decidirá sobre se o 
modelo é adequado ou não, e quais as modificações 
que lhe devem ser introduzidas para que represente 
um equipamento melhor, mais capaz ou mais 
eficiente. O ponto fulcral da modelação é o 
relacionamento entre as ciências fundamentais e a 
especialidade tecnológica. O físico-matemático pode 
desenvolver modelos de alta qualidade, mas não os 
aplica a uma obra tecnológica, nem consegue 
dialogar com os especialistas de um ramo de 
engenharia se não souber a sua linguagem, e não 
compreender os problemas práticos que estes 
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pretendem resolver; o especialista numa área 
tecnológica pode imaginar soluções aparentemente 
possíveis para um problema, mas sem um modelo 
quantitativo não sabe quais as soluções que 
funcionam, como escolher as melhores e como as 
optimizar. O progresso tecnológico baseia-se em 
imaginar soluções novas para problemas práticos e 
encontrar a concretização quantitativamente mais 
eficaz. A modelação é a área de interface entre as 
ciências fundamentais e especialidades tecnológicas, 
e só pode ser abordada eficazmente por aqueles que 
estejam dispostos a atravessar a "fronteira, isto é: (1) 
engenheiros da especialidade que reforcem a sua 
formação física e matemática até atingir um bom 
nível científico; (ii) físicos-matemáticos que estejam 
dispostos a dedicar um esforço substancial ao 
conhecimento e familiarização com a especialidade 
tecnológica a que desejam aplicar os seus 
conhecimentos. O ensino da modelação deverá ser 
feito por especialistas na área de aplicação, quer a 
sua carreira tenha começado pela ciência ou pela 
engenharia; com efeito, só um especialista na área 
tecnológica consegue motivar os alunos de 
engenharia, mostrando a relevância dos aspectos 
físicos e matemáticos para a descrição adequada e 
quantitativa de um processo ou equipamento. Os 
alunos de engenharia interrogam-se frequentemente, 
em especial nos primeiros anos: para que serve esta 
matemática ou física? O ensino da modelação mostra, 
pela experiência, que há da parte dos alunos a atitude 
inversa: nós não sabemos a matemática e física para 
estudar este processo! Muitas vezes sabem-na, mas 
não a sabem aplicar; outras vezes, a apresentação de 
um problema tecnológico importante é uma forma de 
despertar o interesse do aluno para um ramo da física 
ou matemática que ainda não conhecem, motivando- 
os assim para o seu estudo, 


$ 4 - MATÉRIAS DA ÁREA CIENTÍFICA 


Em cada especialidade tecnológica, é possível 
fazer uma lista dos temas de matemática e física mais 
relevantes, que se procura incluir nas respectivas 
licenciaturas ou partes escolares de cursos de pós- 
graduação. Uma comparação de tópicos matemáticos 
e físicos relevantes para várias especialidades 
tecnológicas revela muitos pontos comuns: as áreas 
tecnológicas distinguem-se essencialmente pelo 
domínio de aplicação e linguagem que lhes está 
associada, mas os fundamentos têm numerosas 
afinidades. Com efeito, existe uma gama limitada de 
fenómenos ou efeitos que o engenheiro pode utilizar 
na concepção de equipamentos, e muitos métodos 
matemáticos encontram aplicação nos domínios mais 
diversos da tecnologia. Daqui resulta que é possível 
fazer uma Lista de Tópicos de matemática e física que 


cobrem a maioria das necessidades fundamentais das 
áreas tecnológicas mais importantes na actualidade. 
Esta lista, que se apresenta no Anexo À, consiste de 
36 capítulos, e para cada um destes é indicado um 
exemplo típico de aplicação tecnológica. Na 
realidade, a gama de exemplos podia ser muito mais 
completa, e é apenas para evitar alongar essa lista 
que se restringiu o número de instanciações. Embora 
a lista de tópicos seja algo longa, é muito menor que o 
número de especialidades de engenharia 
multiplicado pelo número de tópicos que cada uma 
usa, visto que a maioria de tópicos têm interesse para 
várias especialidades tecnológicas. De facto, existem 
problemas análogos em vários domínios da 
engenharia, em geral aparecendo sob linguagens 
diferentes e conduzindo a aplicações diversas, que por 
vezes fazem esquecer que se trata de variantes ou 
repetições do mesmo problema. Como exemplos de 
analogias relevantes em várias especialidades 
tecnológicas damos: (i) os campos electro e 
magnetoestáticos em electrotecnia, o escoamento 
potencial em mecânica dos fluidos, a elasticidade 
plana em resistência dos materiais, a condução de 
calor estacionária, a difusão da massa estacionária e 
o campo gravítico clássico em outras especialidades, 
que são todos descritos por equações de Laplace ou 
Poisson; (11) os circuitos eléctricos e mecânicos têm 
analogias, que se estendem a circuitos hidráulicos e 
ramificações acústicas, e têm analogias com cadeias 
de desintegração e reacções químicas; (iii) as ondas 
electromagnéticas, as vibrações de cordas, barras, 
placas e corpos tri-dimensionais, as ondas 
hidráulicas, acústicas e gravíticas constituem outro 
conjunto de problemas análogos. Embora a lista de 
tópicos relevantes em anexo cubra a maioria das 
necessidades das várias especialidades tecnológicas, 
os exemplos dados servem apenas de indicação de 
algumas das muitas aplicações importantes. 


$ 5- INSERÇÃO NOS CURSOS DO L.S.T. 


A lista de tópicos de modelação relevante para 
a totalidade das especialidades tecnológicas consiste 
de 36 capítulos indicados no Anexo A; o grau de 
desenvolvimento apropriado a cada capítulo 
corresponde a cerca de 150 páginas de folhas! A4 
dactilografadas, o que implica que é possível 
leccionar 4 capítulos numa disciplina semestral de 
licenciatura, ou numa disciplina mais intensiva de 
pós-graduação. Portanto, os 36 capítulos de 
modelação correspondem a 9 disciplinas, e seria 
pouco viável incluir a sua totalidade nos cursos 
existentes. É de salientar que, do ponto de vista de 
um certo domínio tecnológico, certos tópicos são mais 
importantes ou têm utilização mais frequente que 
outros, pelo que é razoável dividir as 9 disciplinas em 
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três grupos de prioridades: (i) as disciplinas 
obrigatórias, que versam os temas de modelação de 
uso quase constante, e indispensáveis para a 
especialidade; (ii) as disciplinas opcionais, que 
versam temas importantes para uma boa formação 
do aluno, mas não podem ser incluidas 
taxativamente no plano de estudos, por razões de 
espaço ou outras; (iii) as disciplinas omissas, por 
corresponderem mais a uma formação 
interdisciplinar do que a uma especialização no 
domínio tecnológico em causa. À divisão das 
disciplinas de modelação em fixas, optativas e 
omissas feita acima, tem em vista os cursos de 
licenciatura. No caso de uma pós-graduação, por 
exemplo, a parte escolar de um curso de mestrado, é 
natural que: (ii) as disciplinas optativas na 
licenciatura possam passar a obrigatórias no 
mestrado, se se pretender reforçar a formação dos 
alunos em tecnologias associadas à especialidade: 
(iii) as disciplinas omissas na licenciatura seriam 
incluidas como optativas num mestrado que lhe serve 
de sequência, se este visasse dar uma formação 
interdisciplinar. A classificação das disciplinas em 
obrigatórias, optativas e omissas depende, 
evidentemente, do domínio tecnológico em causa, À 
Tabela de Prioridades para cada Domínio 
Tecnológico, apresentada no Anexo B indica, para 
cada licenciatura e ramo leccionado no I.S.T., quais 
as disciplinas de modelação que, na opinião do autor, 
deveriam ser consideradas obrigatórias e optativas 
ao nível da licenciatura e mestrado. À tabela reflecte 
um esforço para reduzir o número de disciplinas 
obrigatórias, tendo em conta as pressões existentes 
sobre os elencos de disciplinas; não seria dificil, na 
maioria dos casos, justificar a passagem de 
disciplinas optativas a obrigatórias, mas não foi essa 
a opção seguida: tomou-se a via 'minimalista' de 
indicar os tópicos de modelação essenciais a cada 
nível, e apenas esse. À escolha das disciplinas a 
incluir em cada licenciatura ou mestrado e os seus 
programas, poderiam e deveriam ser ajustadas de 
acordo com os docentes da especialidade em causa. 


$ 6 - RESULTADOS DA ACTIVIDADE EM 
CURSO 


As duas tabelas anexas podem ser usadas como 
se segue: (i) do ponto de vista de um certo domínio 
tecnológico, o Anexo B indica as disciplinas 
obrigatórias e optativas, e o Anexo À resume os 
tópicos incluídos em cada uma delas; (ii) do ponto de 
vista de uma formação completa na área cientifica de 
modelação, o Anexo À indica a sequência de 
disciplinas, e o Anexo B os seus principais domínios 
de aplicação. O planeamento da modelação nos 
cursos de licenciatura e mestrado leccionados no 


I.S.T., é bastante menos 'revolucionário' do que pode 
parecer à primeira vista: (i) supõe a manutenção das 
disciplinas do 'tronco comum" de matemática e física, 
que oferecem uma formação de base preliminar à 
área de modelação; (ii) não exclui a existência de 
licenciaturas em matemática e física, às quais a 
modelação deve dar um contributo importante, visto 
que é a área priveligiada de relação com a 
engenharia e tecnologia; (iii) defende que as 
disciplinas de modelação sejam leccionadas por 
docentes da especialidade tecnológica em causa para 
melhor motivar os alunos; (iv) corresponde, em parte, 
ao que já se pratica no 1.S.T., visto que a modelação é 
de facto ensinada, de forma mais ou menos dispersa, 
em todas as licenciaturas e mestrados. O 
aparecimento da modelação como área científica 
distinta verificou-se no Departamento de Engenharia 
Mecânica, traduzindo-se nas seguintes disciplinas: 
(1/ii) a Mecânica Aplicada I e II, respectivamente no 
10 e 20 semestres do segundo ano, que correspondem 
às disciplinas IV e V (do Anexo A); (iii/iv) a 
Matemática Aplicada e Análise Numérica, no 10 
semestre de, respectivamente, o terceiro e quarto 
anos, que correspondem aproximadamente às 
disciplinas 1 e III (do Anexo A); (v/vi) mecânica 
varicional e complementos de matemática 
respectivamente do 10 e 20 semestres do quarto ano 
do ramo de Sistemas, e várias disciplinas de 
mestrado, que correspondem às disciplinas Il e VI (do 
Anexo A). Em resumo, foi acumulada ao longo dos 
últimos dez anos, uma experiência de mais de 6000 
alunos x disciplina, no leccionamento de modelação a 
alunos da licenciatura a partir do 2º ano, e de vários 
mestrados. A experiência adquirida demonstra que: 
(1) os alunos não põem em dúvida a utilidade dos 
fundamentos de matemática e física, se a sua 
exposição for acompanhada de aplicações detalhadas 
à sua especialidade de engenharia; (ii) o ensino de 
modelação pode ser escalonado ao longo do curso de 
licenciatura, dos primeiros até aos últimos anos, 
permitindo atingir um nível de conhecimentos mais 
elevado; (iii) o desenvolvimento da modelação em 
paralelo com as outras áreas científicas da 
licenciatura, evita que os alunos esqueçam os 
fundamentos quando chegam às aplicações, e 
contribui para o bom funcionamento das disciplinas 
da especialidade. 


$ 7 - CONCLUSÕES 

Do que foi exposto, em termos da experiência 
num departamento, ressalta que a área de modelação 
é "horizontal", no sentido que se cruza com as 
restantes áreas 'verticais' da especialização, além 
disso, a modelação conduz à sua própria 
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especialidade, com incidência especial nos fenómenos 
interdisciplinares, de importância crescente na 
tecnologia moderna, e dificeis de abordar por outra 
forma numa estrutura de ensino tradicional. O 
aparecimento de áreas científicas de modelação nas 
várias licenciaturas e mestrados do 1.5.T., no âmbito 
do plano-quadro exposto no anexo B, teria as 
seguintes vantagens: (i) ao nível de cada 
Departamento, encarregado de leccionar um domínio 
tecnológico, permitiria desenvolver os aspectos de 
modelação, essenciais na tecnologia moderna, e 
estabelecer uma relação mais forte com os 
Departamentos de base (Matemática e Física); (ii) a 
existência de uma área de modelação nas várias 
licenciaturas facilitaria o ensino e investigação em 
domínios interdisciplinares, que são as principais 
áreas de progresso tecnológico na actualidade, 
minorando os efeitos de uma divisão demasiado 
rígida, e algo desactualizada, da engenharia em 
especialidades estanques; (iii) o facto de haver um 


elenco de disciplinas comuns a todos os cursos, 
ajustadas a cada especialidade e leccionadas pelo 
Departamento respectivo, criaria um padrão de 
qualidade de ensino uniforme no I.S.T., sem 
contrariar a autonomia dos Departamentos. A 
modelação poderia ser integrada nos vários cursos do 
[.S.T., tomando como ponto de partida a adaptação de 
um conjunto de livros que, apesar de terem sido 
utilizados até agora apenas em engenharia mecânica, 
foram escritos numa concepção interdisciplinar, 
aplicável às outras especialidades de engenharia. De 
entre esses textos, estão completos e em vias de 
publicação dois livros, a "Análise aplicada e campos 
de força”, patrocinado e em fase de impressão pela 
Fundação Calouste Gulbenkian, e "Complex 
functions and potential theory", aceite para 
publicação pela World Scientific, que se incluem na 
bibliografia breve, constante do Anexo C. 


18 TÉCNICA 2 - 88 


ANEXO A 


LISTA DE TÓPICOS DE MODELAÇÃO 
(1- 1X,1-36) 


Consideram-se a álgebra, análise e física geral como 
conhecimentos de base, e indicam-se 36 capítulos, 
organizados em 9 disciplinas, com um exemplo de 
aplicação para cada. 


I- CAMPOS POTÊNCIAIS 


1 - Domínio complexo: O uso de complexos no estudo 
da estabilidade dinâmica de sistemas mecânicos e 
circuitos eléctricos; 

2 - Derivação e integração: uso de funções holomorfas 
para analisar e sintetizar escoamentos potenciais, 
campos electroestáticos, etc.. 

3 - Séries de potências: decomposição de campos 
potenciais em multipolos e o uso de singularidades no 
método das imagens; 

4 - Transformação conforme: escoamento potencial 
em torno de perfis de asa, através de fendas, à saída 
de canais e em jactos livres; 


|I-FENÓMENOS DE TRANSFERÊNCIA 


5 - Funções transcendentes: cálculo do campo 
electroestático de condensadores de placas em várias 
posições relativas; 

6 - Funções generalizadas: funções de influência ou 
de Green para cordas, barras e campos escalares ou 
vectoriais; 

7 - Equações diferenciais ordinárias: oscilações livres 
ou forçadas, com ou sem dissipação, lineares ou não, 
de sistemas mecânicos e circuitos eléctricos com 
vários graus de liberdade; 

8 - Equações às derivadas parciais: fenómenos de 
propagação e difusão em linhas electromagnéticas, 
canais hidráulicos, barras elásticas, tubos acústicos, 
etc. 


II- ANÁLISE ESPECTRAL 


9 - Séries de Fourier: análise de vibrações de motores, 
estatística de processos discretos; 

10 - Integral de Fourier: fenómenos de difusão não 
estacionária de calor, massa ou carpa eléctrica, 
estatística de variáveis aleatórias contínuas; 

11 - transformada de Laplace: estudo da resposta de 
sistemas de controlo de entrada-saída com 
realimentação e acoplamentos vários; 


12 - Métodos aproximados: cálculo aproximado de 
frequências próprias de vibração, modos de 
instabilidade elástica, etc. por via analítica ou 
numérica. 


IV-MECÂNICADOS CORPOS RÍGIDOS 


13 - Álgebra de multiplicidades: equilíbrio estático de 
forças, momentos e tensões em estruturas 
recticuladas e contínuas; 

14 - Transformações de coordenadas: vários grupos de 
transformação de coordenadas e sistemas de 
coordenadas ortogonais adaptados às formas mais 
importantes; 

15 - Corpos e colisões: geometria de massas e 
momentos de inércia de corpos, e colisões envolvendo 
translação e rotação; 

16 - Lei de transformação tensorial: classificação de 
tensores e aplicação ao estudo de trajectórias e 
movimentos giroscópicos; 


V-FÍSICA DOS MEIOS CONTÍNUOS 


17 - Métrica e derivação: aceleração em trajectórias 
curvilíneas e deformação de corpos de formas não 
planas; 

18 - Curvas e superfícies: movimento ao longo de 
curvas e rolamento com escorregamento sobre 
superfícies; 

19 - Gradiente, divergência e rotacional: operadores 
na forma tensorial generalizada, e o seu uso em 
elasticidade com fenómenos térmicos e 
electromagnéticos; 

20 - Operadores de transporte: leis de conservação de 
massa, momento e energia em fluidos dissipativos 
(viscosos, com condução de calor e difusão de massa); 


VI- MÉTODOS VARIACIONAIS 


21 - Extremos e variações: equações de Euler- 
Lagrange para sistemas com N variáveis escalares, 
vectoriais ou tensoriais, dependentes do espaço e 
tempo, e derivadas até qualquer ordem, e com 
constrangimentos; 

22 - Espaço de fases: formulação Hamiltoneana da 
mecânica clássica e dos meios contínuos; 

23 - Geodésicas e curvatura: optimização de 
trajectórias num campo gravítico clássico-relativista; 
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24 - Dimensões e unidades: relações de 
similitude e dimensionais na análise de 
observações e estabelecimento de leis 
experimentais; 


VII- VIBRAÇÕES E OSCILAÇÕES 


25 - Funções auxiliares: funções gama, beta e 
zeta, integrais de Fresnel e fenómenos de 
difracção; 

26 - Equações diferenciais singulares: funções 
de Bessel e oscilações com simetria cilíndrica 
singular; 

27 - Representações integrais: funções de 
Legendre associadas e problemas com 
simetria esférica; 

28 - Operadores de diferintegração: funções 
hipergeométricas e partículas em poços de 
potencial quânticos; 


VIII- ACÚSTICA E ONDAS 


29 - Sistemas de funções ortogonais: 
polinómios de Hermite e transmissão através 
de meios determinísticos e aleatórios; 


30 - Transformadas integrais: 
transformadas de Hankel ou Hilbert, e 
resolução de problemas com fontes; 

31 - Problemas de ordem elevada: vibrações 
acopladas e ondas em meios não- 
homogéneos e em movimento com absorção 
e reflexão; 

32 - Equações diferenciais periódicas: 
ressonância simples e paramétrica em 
sistemas com coeficientes periódicos; 


IX-PROCESSOS 
INTERACTIVOS 


33 - Funções elípticas: solução de 
problemas não-lineares relacionados com 
giroscópios e trajectórias em campos de 
força; 

34 - Fenómenos não-lineares: estudo de 
vibrações de grande amplitude e formação 
e decaimento de ondas de choque; 

35 - Instabilidade e bifurcações: 
determinação de regiões de estabilidade e 
transições entre estados; 

36 - Equações integrais: fenômenos de 
acoplamento sólido-fluido através de 
interfaces, estatísticas de processos de 
colisão. 
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ANEXO B 


PRIORIDADES EM CADA DOMÍNIO TECNOLÓGICO (8 Licenciaturas) 
disciplina indicada na Tabela À 


fenómenos 
transferência 


x 


análise 
espectral 


* disciplina obrigatória no curso de 


licenciatura; 


x disciplina optativa no curso de 
licenciatura, e obrigatória no curso de pós 


graduação; 


IV V VI VII VII IX 
corpos meios métodos vibrações ondas € processos 
rigidos contínuos  variacionais oscilações acústica interactivos 

A) + = E a o 

E x x - - 

= K X x 
: . - x x 

x * = E] x x 

X X x x * o] 

Li + + x X X 

+ & + x m . 

+ * +“ x a 

+ 4 x x E pn 

= * + = É. E 

X K x E = = 

x E x - - - 

x + : + ) E 

X é - - - - 


- disciplina facultativa no curso de 
licenciatura, e optativa no curso de pós- 
graduação; 


A tabela seguinte poderia ser 


revista, tendo em conta os pareceres de 
docentes de cada domínio tecnológico. 
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ABSTRACT 


A matrix formalism for electromagnetic wave 
propagation in multilayered bianisotropic media, 
intended for application either in transmission and 
reflection of plane waves or in guided wave 
propagation, is presented. 

Two specific problems are considered: the 
transmission and reflection of a plane wave when 
propagating across a multilayered bianisotropic 
medium and the derivation of the modal equation for 
an open multilayered bianisotropic waveguide. For 
the first problem, the particular case of an air- 
magnetoplasma interface is studied as an example of 
application of the developed matrix formalism for 
non-guided wave propagation. 


RESUMO 


Desenvolve-se um formalismo matricial para 
meios bianisotrópicos estratificados, para aplicação 
quer a problemas de transmissão e reflexão de ondas, 
que a problemas de propagação guiada. 


ORIGINAL RECEBIDO PARA PUBLICAÇÃO EM 17.02.1988 


Consideram-se dois problemas concretos: a 
reflexão e transmissão de uma onda plana através de 
um meio estratificado bianisotrópico e o 
estabelecimento da equação modal para um guia 
aberto bianisotrópico estratificado. Relativamente ao 
primeiro problema o caso particular de uma interface 
ar-magnetoplasma é estudado como exemplo de 
aplicação do formalismo matricial desenvolvido para 
propagação não guiada. 


|- INTRODUCTION 


Layered anisotropic media have been a subject of 
increasing interest in the last few years due to its 
application either in applied optics or in guided wave 
propagation including millimeter-wave and optical 
integrated circuits, 

The most general approaches to multilayered 
anisotropic media appear to be the 4x4 matrix 
techniques. These techniques have been first 
introduced by Teitler and Henvis [1] and resumed by 
Berreman [2] and Yeh [3]. More recent studies using 
the same basic formulation include the works of 
Schwelb [4], who adopted a transverse resonance 
approach based on impedance matrices for guided 
wave structures, and Krowne [5], whose analysis in 
the Fourier transformed domain is particularly 
suited for multilayered structures such as 
transmission lines and printed antennas. Berreman's 
and Yeh's formalisms have been compared by Lin- 
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Chung and Teitler [6] who developed a unified 
framework for ellipsometric studies. Applications of 
the 4X 4 matrix techniques to ellipsometric and 
diffraction analyses concerning anisotropic media 
have been also recently presented [7, 8]. 
Bianisotropic media have been considered in 
References 2 and 5. 

In this paper, a general matrix formalism 
applicable to electromagnetic wave propagation in 
bianisotropic media is also developed. 

In Section 2 the basic matrix field equations are 
derived. The mathematical treatment is formally 
similar to other published analyses, namely to 
Berreman's and Krowne's quoted works. However, in 
Berreman's approach, only non-guided wave 
propagation was considered as the focus of interest 
was optics, restricting the directions of incident 
waves to a plane perpendicular to the layered 
medium. On the other hand, Krowne's analysis is 
devoted only to open and closed guided wave 
structures. The mathematical treatment developed 
in Section 2 is more general, in the sense that either 
non-guided wave propagation or guided wave 
structures may be considered. For this purpose, all 
three rectangular components of wave vectors are 
considered, in spite of the homogeneity of the 
structure in the horizontal plane (Fig. 1). Hence, for 
ellipsometric analyses any direction of incidence 
may be considered and for guided wave propagation 
along the z-direction one only has to make k,=0 in 
all equations derived in Section 2 (see Fig. 1). The 
analysis presented in Section 2 includes general 
expressions for the electromagnetic fields, based on 
the eigenvalue-problem for the derived 4x4 coupling 
matrix which results, for each layer, from the 
application of Maxwell's equations. 

In Section 3 the problem of reflection and 
transmission of plane waves at an air- 
magnetoplasma interface is solved. This is a simple 
example, which has the merit of illustrating the 
application of the general theory developed in the 
previous Section and providing a comparison with 
available results [9]. 

In Section 4 a generalization of the problem 
studied in Section 3 is carried on by means of the 
definition of a transfer matrix for the multilayered 
medium. 

Finally, in Section 5, the theory presented in 
Section 2 is applied to the derivation of the modal 
equation for open bianisotropic multilayered 
waveguides. Although this modal equation is quite 
general, particular attention is paid to lossless 
crystals at optical frequencies, namely to the 
particular form of its coupling matrix. One should 
note, however, that - as far as the authors are aware - 
this constitutes a new approach to bianisotropic 
layered waveguides. 


Fig. 1 - Layered medium with n - 2 interior layers of 
total thickness t. 


2 - BASIC MATRIX FIELD EQUATIONS 


The stratified medium under consideration may have 
an arbitrary number of anisotropic layers as depicted 
in Fig. 1. 

Electromagnetic fields are assumed to have 
explj(wt-kyy-k;z)] dependence. Since the whole 
layered medium is homogeneous in the yz plane, ky, 
and k, remain the same throughout all strata. 

Wave vectors k have the general form 


k=ht+k,y+k,2 (1) 


For ellipsometric problems, given k, and k, one has 
to calculate the x-components hgfts) for the sth-layer 
with Isssn and where n is the total number of 
layers. The o-subscripts account for all possible 
values of hts) in each layer. For guided wave 
propagation along the z-direction, one has k,=0 due 
to the uniformity of the structure along the y- 
direction and one has to calculate k, and all the 
transverse wave numbers hs), 

The complex Maxwell's equations in frequency 
domain 


VxE=-jwB 
(2) 
VxH=jwoD 


may be written in matrix form for the (x,y,z) frame of 
Fig. 1 as | 
Vx 10 1 | .e=jwf 
(3) 
1 0 
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where e and f field vectors are defined as 
e= IE É f= |D 
(4) 


H B 
and where O and 1 are, respectively, the 3x3 null 
and unit matrices. 

For the sth-layer, the constitutive relations in 
irequency domain are 


D= eq£ (00). E +E().H 


(5) 
B=0"(0).E+ Hop co). H 
where [9, p.68]t for each tensor a in (5) 
a(w)= | ale !“*dt 
í (6) 


In Egs. (5) the labelled tensors are independent of the 
space-coordinates for each layer, i.e. each layer is 
homogeneous. The E and ( tensors are called 
magnetoelectric tensors [9, p.62] and are nonzero in 
optically active media [2, 5]. A medium characterized 
by Eqs. (5) is said to be bianisotropic [10, p.71. 

If one defines a 6X 6 matrix for the (x,y,z) frame 
ofFig.1 as 


E) | 
Ms ae e ade ES 


(8) (s) 
[4 


Hop 


the constitutive relations in Egs. (5) may be 
condensed in the form 


E=M, 
- (8) 
Defining a 3X 3 matrix operator 
R,=| 0 jk, -Jk, 
(9) 
| d 
Jk, O — Fa 
d 
Jk, FR 0 


(i nara superscripts (s) are Ena away 
to simplify the notation where there is no risk of 
confusion. 


Eq. (3) may be written in compact form as 


R.e=jwf (10) 
where 
R=- 0 R, 
— R, O 
(11) 
isa 6x6 matrix operator. 
Using Egs. (8) and (10), one may write 
R.e= joM' e. 
(12) 


Eq. (12) shows that ey and es components of e (i. e., Ey 
and H,) can be algebraically expanded in terms of the 
other field components using rows 1 ard 4 


b 
Jh,es—JR eg=J0 ba Me, 


1=1 (13) 


6 
—Jk,e pk, e =jo 2 M, e. 


= 1 


Rows 2,3,5, and 6 of Eq. (12) lead to the following 
first-order linear differential equations 


de, 6 
—Jk ,e,— ses = jw à My e, 
= 


de, 
det io M, e 


(14) 


de, b. 
Je, + de = jq) > Mae 


1=1 


de, | 
=Jk e;— E = jw > Met. 
t=1 

Eqs. (13) and (14) are similar to Eqgs. (15) and (21) of 
Reference 5. 

In order to handle more easily the boundary 
conditions of the layered medium one introduces the 
field vectors 
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= | €p,€ 


(15) 


W = 


€,,€.,0,4€ ã 
2" 63º 5 


where the symbol "T” stands for transpose. In Egs. 
(15) & includes the normal-to-the-interfaces field 
components whereas 'F includes the tangential-to- 
the-interfaces field components. The special 
arrangement chosen for the 'F field vector will be 
understood later. 

Eqs. (13) can now be written in terms of & and 'F; 


— qts) 7 
p(x) =G P(x) (16) 


where G is a 2x4 matrix whose elements are 
(lI=i<4): 


1 € 
G,= A | MM A DM AM + As] 
(17) 
n 1 
Goj= A MM pt A)-MyM+tA,) 
where 


A=MMa-M,M 
11º dd l4 4i (18) 


and A; are the elements ofa 2X4 matrix given by 


(19) 


In Eqgs. (17) p subscripts are related to i subscripts 
according to Table 1. 


Using Eq. (16) to eliminate e; and eq from Egs. 
(14), one obtains 


d 
pd he JC! xp (x) 


(20) 
where C is a 4X4 coupling matrix whose elements 
are 


| k 
vão = 
Cu=a|M AM Sa (Mat % Je 


k , 
| z 
Cy=o|M,+M,6,+(My+ = e 


k 
Cy= “o|MGAM 409 +(Mg- “e 


RE E 
Cg=-0[M,+M46,+(My E es 


with Isis4; p subscripts are related to i subscripts 
according to Table 1. 

The matrix equation (20) constitutes a linear 
homogeneous differential system with constant 
coefficients for the sth-layer. The expression for 
Wí(s/x) can be found by means of the corresponding 
eigenvalue-problem related to the Cts) coupling 
matrix. 

The eigenvalues of the coupling matrix for the 
sth-layer constitute all the possible x-wavenumbers 
ho!s). These eigenvalues are determined by solving 
the characteristic equation [11, p.200] 


dei RL = 0 | =0 
(22) 
When expanded, this equation reduces to a Booker 


quartic equation [9, p.200Jtii) 
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hº-B h'+Bh!-Bh+B,=0 


(23) 
with four roots hs, i. e. 1=S0<=4. Hence, the 
diagonalization of Cts) is possible if it has four 
linearly independent eigenvectors [11, p.735] ug!s) 
which may be found according to the equation [11, 
p.200] 


ht -C ut = () 
o o 


(24) 
Except for the limit case of an isotropic layer there 
are always four distinct eigenvalues hg!s), thus 
making possible the diagonalization of Cts) [11, 
p.719]. However, even in the isotropic case, where 
each eigenvalue has a degree of multiplicity two, it is 
possible to choose two linearly independent 
eigenvectors for each degenerate eigenvalue, thus 
making always possible the diagonalization of the 
coupling matrix. 

The coefficients B; of the Booker quartic equation 
(23) are conected to the eigenvalues h, by the well- 
known root-coefficient relations [11, p.717] 


B=2h, 


= Ah, 


v>0 


(25) 


> hohh, 


p>v>0 


B,= [1h 


where 0=1,2,3,4 for B and B,0=1,2,3, v=2,3,4 for 
Bb, ando=1,2,v=2,3,p=3,4 for B,. The coefficient Bo 
is also equal to the sum of all the principal minors of 
order o in C [11, p.718]. Therefore, for o=1 one 
obtains the trace of C, i.e. according to the first 
equation in (25) 


(ii) Well known frei studies of wave a is in 
the ionosphere 


When the eigenvalues of the coupling matrix observe 
antysimmetry, i. e. ifone may write 


(27) 


h=-h 

4 2 
the Booker quartic equation reduces to the 
biquadratic form 


(28) 
and B/=By=0. 
For the sth-layer, a diagonal matrix 
(8) ts) pls) (8) p(8)| 
H'º=dia en ho ER EA 
(29) 


may be obtained using the standard expression 


H'º — | uts) | — Eqts) U'” 


(30) 
where Uts) is the diagonalizer matrix for C(s),i.e. 


+ (8) (5) (Cs) (s) 
U = [uu ul? ul? 
(31) 
according to Eq. (24). 
Introducing a new vector 
E (7 (8) (8), (8) (8) T 
a (w)=|a, (x),a, (x),a9 (x),a, | 
(32) 
related to the w(s) (x) field vector by 
(33) 
then Eq. (20) reduces to 
d 
Es a! (x) sumi - IH a! (x) 
(34) 
whose solution may be written in the form (1=0<4) 
18) 
a = 418) a SER 
O o 


(35) 
where Ag's) is an amplitude referenced to the 
entrance coordinate of the sth-layer, x=x,.1 (see 
Fig.1). Therefore, the w's) (x) field vector for the sth- 
layer is 
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é —jh 
wa a? Al (8), o s-—l 


(36) 

The application of the boundary conditions is 
now a straightforward task, since it reduces to 
imposing the continuity of wy on the several 
interfaces of the multilayered structure. 

The formalism of this Section is a generalization 
for bianisotropic media of a previous work of the 
authors[12]. 


3 - REFLECTION AND TRANSMISSION OF 
PLANE WAVES- AN EXAMPLE 


In this Section the matrix formalism developed 
in the former Section is applied to an example of 
reflection and transmission of plane waves at an air- 
magnetoplasma interface as depicted in Fig. 2 [9, 
pp.287-291]|. Medium 1 is the air and medium 2 is a 
magnetoplasma. The interface is located atx=0. 


X 


a 
* 
* | 


Ci 


% 
Cd 


Fig. 2 - Boundary - value problem of an air 
magnetoplasma interface. The plane of incidence is 
the xy plane (k, = 0). The externally applied 
constant magnetic field is z - directed and the wave is 
propagating normal to the applied field. 


The externally applied constant magnetic field is z- 
directed, i. e. B.=B 2. Hence, for s=2 (medium 2), 
Egs. (5) reduce to 


D= ee). E 
(37) 


B=yH 


where [9, pp.264-5] 


ET, Je, Õ 
(38) 
Je, E, Õ 
É O o 
and 
o 
Ó 
e (o)=1— ; 
1 a 
4979 
e (m)= — : 
l-E, 
(39) 
8 (w)=1 milho 


with the frequency dependent coefficients X, and Y, 
given by 


Er ta Up 
Em od ae 
(40) 


and where wp is the plasma frequency and wp the 
gyrofrequency of the electrons. When Bç=0, 
Wp= Y,=2,=0,, 


(41) 


and the plasma behaves like an isotropic medium. 

The plane of incidence is the xy plane which is 
perpendicular to the air-magnetoplasma interface, 
thus making k,=0. According to Fig. 2, one has, with 
ko = w(£oHo) 


k =k sinb. 
y É) t 
(42) 
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where O is the angle between the incident wave 
vector k; and the positive x-axis. 


In medium 1, wy'!) may be decomposed into WU l)=A 0 RE 1 da 
incident wy; and reflected wy, feld vectors. According À | 
to Egs. (21) one has for the air region 
CPE 4 , Do 
l ; te 
(43) h, 
0 C | 
' (48) 
where 
1 0 
h2 
C=|0 =] €=] é =p E 
1 DE r ( da O 
' (44) OH, 
2 
ho and 
(DE Ú — O 
UU OH, 
and with 
jh x Jh x 
a ã 
Wp,(0)= B, 0 e +B, 1 e 
hê=hº- po. WE 
a ) y = 0 | A 
(45) a 
For this isotropic region there are only two distinct a 9 
eigenvalues for C"': tha. From Egs. (42, 45) one has (49) 
h =k cos0 1 0 
a Ú E a 
(46) 
For the incident (reflected) field vector w, (wW,) one h. 
chooses h!!!= + ha(-ha), Hence, the incident and the | 0 
| SMA | wH 
reflected wave vectors are given, respectively, by o | 
where A;,, B, define the perpendicular polarization 
k=hX+ky and Ay, Bythe parallel polarization. 
? For the magnetoplasma, according to Egs. (37- 
(47) 38), one obtains from Egs. (21) 
k=-h%+kY C“=| C 0 
r a y e 


(50) 
An appropriate choice for Cl eigenvectors(i!) gives 
(see Eq. (36)) 
(11) There must be four eigenvectors forming an U 
orthogonal set (sec Section 2). 
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where 


e e hi ho 
pi, A 
e, WE E, 
(51) 
E Eh e 
—— |eó —& Te b 
lej=) qe 
Bicas E 
and 
= 0 —wN 
(52) 
his ki 
À 0 
tw 


In this gyrotropic case the four roots of the 
Booker quartic equation (23) are th, and the, with 


(53) 
hº= [2 E) 2 p? 
x o Y 


therefore observing antisymmetry (see Eq. (27)) and 
making null the trace of the corresponding coupling 
matrix, i.e. 


(54) 


a result which agrees with Egs. (50-52), One should 
note that he (ho) is associated with the extraordinary 
(ordinary) wave as it does (not) depend on B.. Taking 
into account that for medium 2 only the eigenvectors 
corresponding to eigenvalues of positive value can 
propagate into the magnetoplasma, the transmitted 
field vector is 


-=Ih x 
| pe * 
L 
e 
(55) 
O 
O 
where 
| Poa 
ho oe 2 
aa Sra 
9 wu E e he k, 
(56) 


according to Eg. (24). Hence, in the magnetoplasma, 
there are two transmitted wave vectors 


k=hx+ky 
o o y 
(57) 
a A A 
k,=h,X+ky 


with ho, he>0 as depicted in Fig. 2(iv). 
The boundary conditions in the x=0 plane reduce 


LO 
wp (0) +wy (0)=wy (0) 
(58) 
Therefore, 
C, =|1, 0 JA, 
(59) 
C, 0 uy Aj 
and 
B = E O A 
(60) 


B b Ita 
with, for 6, =n/2, 
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(61) 


Pd E 
| I4+Tvh/wr 
e UU E] 


Hence, according to Eqgs. (59-60), the perpendicular 
polarization excites the ordinary wave while the 
parallel polarization excites the extraordinary wave, 
existing no coupling between the two polarizations in 
medium 1. 

When B,=0, Eq. (41) is valid. In this case, 
according to Egs. (53) 


2 2 2 z 
h'=h'=ek —k 
7 e 1) y 


(62) 
and the second equation of (61) reduces Lo 
E 2 
TO 1+eh /h 
io 
(63) 


In the limit case ofe-—1 one gets the expected results: 


Ee | | À L, =4=1 
(64) 
For the particular case of normal incidence 


(0;= 0) one obtains from Egs. (61) 


1/2 


(65) 

However, Eqgs. (61) cannot be applied for grazing 
incidence (0,=n/2) as h,=0. Nevertheless, a 
different analysis is easily carried on in this case: the 
ordinary (extraordinary) wave cannot be excited as 
the continuity for H, (E,) on x=0 is not observed. 
Therefore, v; =1/=0. Moreover, due to the 
uniformity along the y-direction, "| =Tj=0. 


(iv) One should bear in mind the difference between 
the ordinary wave vector ko and the vacuum wave 
number k,;=cw(E pu 


4. TRANSFER MATRIX 


For the layered medium of Fig. 1 one can define a 
transfer matrix relating wtn) (t) to yí(0). 
The boundary conditions to be considered are 


ts +1) (x) 
E) 


py (1) =1wy 
(66) 


with |I<s<n-l. However, according to Eg. (33) these 
boundary conditions may be also written as 


5 s+1,s 8 
(67) 
also withl<s<n-1 and where matrices B; (2<i<n, 


I=j<n-l)are given by 


B = 
t) 


a oa 


(68) 
One can introduce transmission matrices 4'!s) 
with 2<s=<n-l so that 


(8) 418] is) 
a u)=0 .a (x) 
(69) 
where 's' are diagonal matrices given by 
ii a ELA -jhg 8a 
p E diag É ' Ê i 
(70) 
DO Veg é 0082 
d 4 
e e 


and where Axts) is the thickness of the sth-layer, 1. e. 


AX =x — X 
5 s—l 
(71) 


with 2=s=n-l., 


From Eqs. (67, 69) one has 


tn) (1) 
a ()=Qa (0) 
Q (72) 
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where 


Applying Eqg. (33) for s=1,n to Eg. (72) one 
finally obtains the transfer matrix T defined by 


yo (g=7y'D (0) 
(74) 
as 


u'? - | 
(75) 


One should note that IP) (0) is the total field 
vector and not the incident field vector "P;(0). In fact, 
if, (0) is the reflected part of WC) (0), Eq. (74) leads 
to 


vinho qin-l) (2) 
T=U""2 9" "2..B460º. 


wp ()= Ty, (0) + wy, (0)) 
(76) 

where "FP, (t) = Wtn)(t) is the transmitted field vector. 

For the particular case of n=2 as in Fig. 2 
(Section 3), the matrix T reduces to the unit matrix 
and thus Eq. (76) reduces to Eg. (58) in Section 3. For 
a slab between media 1 and 3,1. e. n=3, the transfer 
matrix Is 


T=u2 q. u'? — 1 


(77) 


5. GUIDED WAVE STRUCTURES 


In this Section guided waves along the 2- 
direction in the open waveguide of Fig. 1 will be 
considered. Due to the uniformity along the y- 
direction k,=0 and k, is the longitudinal wave 
number to be calculated. 

In the first part of this Section some 
characteristics concerning the matrix equations (16, 
20) at optical frequencies will be analysed due to the 
practical importance of optical waveguides. In the 
second part, the modal equation for open 
multilayered bianisotropic waveguides will be 
derived. 

IH one disregards optical activity, E(s!= 0s)=0 [2, 
51. For optical frequencies one may also write uís)= 1 
(13, p.132]. Hence, the matrix M(s) defined in Eg. (7) 
reduces to 


(78) 


Moreover, if only lossless (uniaxial or biaxial) 
crystals are considered, the dielectric tensor is 
symmetric [9, p.182]. Therefore 


(79) 


where, as usual in this paper, (s) superscripts are 
dropped away to simplify the notation. 
Accordingly one has from Egs. (17,21) 


E” E k, 
G=|-— 0 -—-— 
E WE 
xx xx [i | 
(80) 
k, 
o ns O O 0 
won, 
and 
C= 0 wu, 0 0 
A ko ea, Na Xy 
O -— A k 
wH E cw " 
0 XI LX XX 
A” - à eus ; hi-e kh 
E. * E * qpE 
II II ÚO LI 
we WE, E, 
-—=S A OO =— À — — k 
E JE E xy E 2 
LX XX J 
(81) 


where Ajj is the algebraic complement of the e;; 
element of matrix £ given in Eq. (79). 
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For the equatorial configuration [4, 12] in which 
Exy=Ey2=0, Ay,=0 and the coupling matrix may be 
written in the form 


(82) 


thus making possible the existence of TE and TM 
modes, provided that all the layers of the waveguide 
belong to the equatorial configuration, i.e. 
exy9=eyils!=0 for every s (|ss<sn). (Y) Then, the IP 
vector may be decomposed into two field vectors 


IE, H i 
Z y 
(83) 
therefore accomplishing the separation of Eq. (20) 
into two independent equations - one for TE modes 
and another for TM modes. Similarly one has (see 
Egs. (15)) 


rw 


TM” 


Pre=H,» Oyy É, 
(84) 

One should note that the chosen arrangement for 
the PY vector defined in Egs. (15) is finally 
understood if one bears in mind Egs. (83). 

For the particular cases of longitudinal and polar 
configurations |4, 12], Eqgs. (80-81) may be further 
simplified. In fact, for these two anisotropic 
configurations, the eigenvalues of the coupling 
matrices observe antisymmetry (see Eqgs. (27)), 
therefore making null the trace of these matrices 
(see Eq. (26)): 


fr 1C]j-0 
(85) 
For the unaxial longitudinal configuration, 


Ex2=Ey,=0, Etr, Ag=epes (vi) and the 
transverse wavenumbers are given by |12|] 
2 2 p? 2 4 : a 
hê=e bé ko, h=>(eh—k? 
o Log e e pu 2 
XX 
(86) 


For the unaxial polar configuration, txz=Exy=o0, 
Eyxx= EL, Ôgx=ep es (VM!) and the transverse 
wavenumbers are given by [12] 


(87) 
Next, let us find the modal equation for the open 
waveguide on Fig. 1. For the general case, k, is 
complex as well as the transverse wavenumbers. 
Henceforward, it will be assumed that the 
eigenvalues for the exterior layers (i. e. for s=1,n) 
are such that 


Rh )>0 Reh)>0 
(88) 


Reh)<O Reh )<O 


Therefore vectors all,,, and at?) (;, introduced in Eq. 
(32) may be decomposed into the two following parts 
(s=1,n): 


es“ to=|2 0 (8 to" 
(89) 
where 
is] 
— Jh ix = ) 
(s) — als) 1 s-—1l 
a, ()=A, e 
(90) 
(8) 
-Jho (x-20 1) 
e 
and 
(8) 
Jha tx- 
a Ls) (x)=A (8) e J g = | 
(91) 
(s) 
— jh (x—x ) 
e 4 8 — | 


A 's/ and Aístare 2X2 diagonal matrices whose 
elements are complex amplitudes. While a ts) (x) 
represents waves propagating in the positive x- 
direction, as! (x) represents waves propagating in 
the negative x-direction. For guided wave 
propagation in the open multilayered waveguides of 
Fig. 1, one must have 
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ig ()=a!P (x)=0 
(92) 
Therefore, Eq. (72) in Section 4 may now be written 


as 
as 01] =19% Vl 0. 
(93) 
(1) 
1) Q,, Vo a (0) 


where Qi are 2x2 submatrices of Q given in Eq. 
(73). A nontrivial solution for Eg. (93) is found if, 
and only if 


det |IQ.]=0 
Voz (94) 


which is the searched modal equation for the open 
waveguide of Fig. 1, allowing the calculation of the 
longitudinal wavenumber k,, which is, for the 
general case, a complex number. 

In spite of the fact that the derivation of Eq. (94) 
was made bearing in mind waveguides used in 
integrated opties, the method is quite general in the 
sense that it is applicable to open multilayered 
bianisotropic waveguides. 


6. CONCLUSIONS 


In this paper a matrix formalism suitable for the 
analysis, in frequency domain, of multilayered 
bianisotropic media is developed. The feature which 
distinguishes this analysis from other similar 
methods is its flexibility in handling non-guided 
wave propagation as well as guided wave 
propagation in multilayered bianisotropic 
structures. 


As an example of application of the developed matrix 
formalism, the boundary-value problem of an air 
magnetoplasma interface was solved. A 
generalization of this boundary-value problem for a 
medium with an arbitrary number of bianisotropic 
layers was also presented, 


(vi) For the uniaxial crystal, edge) is the ordinary 
(extraordinary) relative dielectric constant 


Finally, some relevant characteristics of guided 
wave structures at optical frequencies were analysed 
and the modal equation for an open multilayered 
bianisotropic waveguide was derived. Although only 
open waveguides were considered, the matrix 
formalism presented in this paper can be also 
applicable to closed waveguides. However, this 
special case will be presented elsewhere. 

The method presented in this paper covers a wide 
range of phenomena related to electromagnetic wave 
propagation in time-dispersive bianisotropic media 
and is particularly suited for computer calculations 
concerning multilayered structures. 
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